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Abstract 

In practical combustion systems, the assumption of chemical equilibrium 
is not always appropriate. This is especially true in regard to predicting 
various critical combustion processes including ignition, flame stability, 
smoke and gaseous emissions. In order to estimate combustor performance 
accurately, it is important to calculate combustor internal profiles of 
temperature, species and velocities to reasonable accuracy. This, in turn, 
requires realistic assumptions concerning the turbulent transports of heat, 
mass and momentum and their interaction with chemical reactions. 

Conventional combustion models assume the scalar transport process can 
be described adequately by an assumed two/three parameter probability density 
function (pdf) and that the transport of heat and mass is essentially identical. 
These models give reasonably well correlated results for simple flows. However, 
their applications to complex flows involving strong recirculation and swirl 
do not appear promising. Since such flows are of common occurrence in any 
combustion systems, straight forward’ extension of conventional combustion 
models to combustor flow calculations is not adequate. 

The present approach is formulated to investigate some of the problems 
that are common to both reactive and non-reactive variable density flows. 

Such problems are mass transport in a variable density flow, validity of 
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constant density turbulence models, relative merits of Reynolds versus Favre 
average and calculated versus assumed pdf for the scalar. The approach solves 
for both Reynolds and Favre averaged quantities and calculates the scalar pdf. 
Turbulent models used to close the governing equations are formulated to 
account for complex mixing and variable density effects. In addition, turbulent 
mass diffusivities are not assumed to be in constant proportion to turbulent 
momentum diffusivities. 

The governing equations are solved by a combination of finite-difference 
technique and Monte-Carlo simulation. Some preliminary results on simple 
variable density shear flows are presented. The differences between these 
results and those obtained using conventional models are discussed. 
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MOTIVATION 


• UNDERSTANDING OF TURBULENT 
REACTIVE FLOWS 

• MODELING OF TURBULENT REACTIVE 
FLOWS 

• COMBUSTOR FLOW MODELING 

• DEVELOPMENT OF SUITABLE COMBUSTOR 
DESIGN GUIDES 


TURBULENT REACTIVE FLOWS 



FLOW (VELOCITY, SPECIES 





MODELS) 
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MAIN PROBLEM IN COMBUSTION MODELING 


NON-REACTIVE PLOW MODELING 

• CONTINUITY 

• MOMENTUM (TURBULENT MOMENTUM FLUXES, - pu'.u'. 

^ J 

• ENTHALPY (TURBULENT HEAT FLUXES, , ETC.) 

• SCALAR (TURBULENT SCALAR FLUXES, , ETC.) 

• FLUCTUATING MASS FLUXES (p'u' , ETC.) 

• QUESTION OF REYNOLDS VS. FAVRE AVERAGE 

• VALIDITY OF CONSTANT DENSITY MODEL FOR 
VARIABLE DENSITY FLOWS. 

REACTIVE FLOW MODELING 

• EVALUATION OF MEAN FORMATION RATES 

- R ($, h, p, ) 1 R ($, h, p, ) 

- FAST/FINITE CHEMISTRY 

• VALIDITY OF CONSTANT DENSITY TURBULENCE MODEL 
FOR REACTIVE FLOWS 


, ETC.) 
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ASSUMPTIONS 


• FAST CHEMISTRY 

(IF IT MIXES, IT REACTS) 

• ONE STEP FORWARD REACTION 

(FUEL + OXIDANT - PRODUCTS) 

• TURBULENT TRANSPORT OF DIFFERENT 
SPECIES IS THE SAME 

(TURBULENT SCHMIDT NUMBER 
FOR ALL SPECIES IS IDENTICAL) 

• EDDY DIFFUSIVITIES FOR MASS AND 
HEAT ARE IDENTICAL 

(LEWIS NUMBER EQUALS UNITY) 


CONSERVED SCALAR, f 

f = MIXTURE FRACTION 

= MASS OF SPECIES 1 
MIXTURE MASS 

RELATED TO OTHER CONSERVED SCALAR, 
SUCH AS ENTHALPY, h, BY 

h = (h fu “ + h ox 
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EQUATIONS SOLVED IN CONVENTIONAL REACTIVE 
FLOW MODELING 

. MEAN CONTINUITY 

• MEAN MOMENTUM (CLOSURE VIA TURBULENCE 

MODELING) 

• MEAN CONSERVED SCALAR (CLOSURE VIA 

TURBULENCE MODELING AND SCHMIDT NO.) 

. MEAN VARIANCE OF CONSERVED SCALAR 
(F 7 , WHERE f' = f-f) 


CONVENTIONAL APPROACH 


• SOLVE EQUATIONS USING FINITE 
DIFFERENCE SCHEME 


• ASSUME A TWO-PARAMETER PDF 
FOR f TO ACCOUNT FOR SCALAR 
FLUCTUATIONS 


• DETERMINE THE PARAMETERS BY 
REQUIRING FIRST AND SECOND 
MOMENTS TO AGREE WITH I AND 
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PRESENT MEASUREMENT AND PREDICTION CAPABILITIES FOR TURBULENT REACTIVE FLOWS 




Measurement Capabilities 

Prediction Capabilities 



Physical 

Quantities 

Probe 

Data 

Optical 

Data 

Difficult 

? 

Assumed 

Pdf 

pdf 

Transport 

Difficult 

? 

Remarks 

Flow 

Field 


X 



X 





X 

X 


X 



Hot-wlre/LDA 

technique 


mm 

X 

m 

X 

v,w u 




Hot-wlre/LDA 
technique;Need 
to Improve model 


mm 




a 


X 

« 


h'uj 

X 



X 


X 

II 


Pdf of 

X 


X 

(Major 

Species) 

X 
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Raman scattering 
tech, need to 
improve model 

Combustion 

Field 

*a' U i 
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X 
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(Model) 
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X 
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meas,; Rayleigh 
scattering 
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CONVENTIONAL APPROACH TO 
COMBUSTION MODELING 


• REDUCE PROBLEM TO EQUIVALENT 
SCALAR MIXING PROBLEM 











































LIMITATIONS OF CONVENTIONAL APPROACH 


. INFINITELY FAST CHEMISTRY 

. ONE/TWO STEP REACTIONS ONLY 

. SCALAR FLUCTUATIONS DETERMINED BY 
ASSUMED PDF 

. ALL SPECIES DIFFUSE AT SAME RATE 
• HEAT/MASS DIFFUSE AT SAME RATE 
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PRESENT APPROACH ADDRESSES THE 


FOLLOWING PROBLEMS 

• PDF OF SCALAR 

• MASS DIFFUSION IN TURBULENT FLOWS 

• REYNOLDS OR FAVRE AVERAGE 

• CALCULATION OF FLUCTUATING 
MASS FLUXES 


NOMENCLATURE 

U^ , $ , etc. 

U.^ , etc. 

Uj , 0, etc. 

* 0*» P’ » etc. 

F <0 ) 
a 

< > 


FAVRE AVERAGE 
REYNOLDS AVERAGE 

FLUCTUATING PART OF FAVRE DECOMPOSITION 

FLUCTUATING PART OF REYNOLDS 
DECOMPOSITION 

PDF. OF 0 
TIME AVERAGE 
INSTANTANEOUS QUANTITIES 
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GOVERNING EQUATIONS 


• CONTINUITY 


d 

/ — « 


P F 

dx . 

\ 

i 


MOMENTUM 

d 

r - „ 


p u . 

dx . 

L 

l 


PDF OF 

Mai 


3P 

dx 


U. 

l 



1 d 

f r 9F(? a ) ' 

1 - — [ 

P 0X. 
^ 1 

L dx. 

l 

1 30 L 

a 


F (0 ) S(0 ) 


+ E (0 , x.) 

a 1 
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TURBULENT CLOSURE 


TWO-EQUATION PLUS NON-EQUILIBRIUM ALGEBRAIC STRESS MODEL 


P u. 


at a r_ k _ . . ak i 

= |C _ p <u, u.> | 

3x. 3x. L e J dx. J 


“ P <n, u > iHi 
J 3x. 


- pc 


pc + <£^i> £L 

P 3x . 
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MODEL FOR p AND <p'u'> 


CONSIDER MIXING OF TWO FLUIDS WITH DIFFERENT DENSITIES. 


P 

U. = U. + a. i> <u. 0 > 
1 i 2 1 


<p'u*> = - a. p <u.0> 

l 2 l 
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$ FROM PDF F(0 ) AND <u.0> FROM TURBULENCE MODEL 


MODEL FOR <u^0> 


NON-EQUILIBRIUM ALGEBRAIC FLUX MODEL 


3$ 3U^ p 30 

- <u £ u^> - <nj 0> + < _ > r 0 


dx . 
J 


dx 


j 


p dx. 
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SOLUTION PROCEDURE 


• FINITE DIFFERENCE SOLUTION OF ALL GOVERNING EQUATIONS 


• OPERATOR-SPLITTING TECHNIQUE APPLIED TO PDF EQUATION 


ASSUME F(0 ) = F*(P ) + F'(0 ) 

a a a 

F* AND F' GIVEN BY SOLVING 


0F* = 1 _0_ r 3F* I 

* 3x. p dx, L dx J 

i r i i 


dF* 

dt 


E = 


2° H t J F* 
- 2f t 



(0' q ) p* (20 

F* (0 ) 
a 


+ E(0 , x.) 

a i 


- 0' q ) d0' 


a 
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